Classical stability of supercurrent in one dimension: a numerical study 
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We report results of a classical simulation of thermal phase slips, and the associated relaxation 
of supercurrent, in a ring-shaped one-dimensional superfluid. We find that the classical relaxation 
rate vanishes in the uniform limit. This leaves the quantum relaxation, with momentum transfer 
to phonons, the only mechanism of supercurrent decay in the uniform system. In the presence of 
a smooth periodic potential, classical decay becomes possible, and we identify a family of moving 
critical droplets that can mediate it. 
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I. INTRODUCTION 

Decay of supercurrent in one-dimensional (ID) super- 
fluids and superconductors is a problem of longstanding 
interest. In 1967, Little ,lj discussed the role of "large" 
thermal fluctuations in thin superconducting wires and 
observed that if, as a result of such a fluctuation, the or- 
der parameter $ vanishes at some point, the phase of ^ 
can unwind, leading to a nonzero resistance. Such pro- 
cesses became known as thermally- activated phase slips 
(TAPS). Using thermodynamic arguments, Little had 
pointed out that at low temperatures the rate of these 
processes is exponentially small. 

A more detailed theory of TAPS, based on the Ginz- 
burg-Landau (GL) equations, was constructed in subse- 
quent works [2| and became known as the LAMH the- 
ory. At first sight, it may appear that the LAMH theory 
is essentially an application of the conventional droplet 
model to the specific case of a superconducting wire. In 
this approach, one identifies a critical droplet as a sad- 
dle point solution to the classical equations of motion (in 
our case, the GL equation). Dynamics along the nega- 
tive mode of this saddle point is supposed to mediate a 
transition between states that differ by one unit of the 
winding number. 

On a closer inspection, however, one notices that the 
order parameter does not vanish on the LAMH saddle 
point, except for initial states with exactly zero current. 
This is because the only current present in the GL theory 
is the London supercurrent / s , which is proportional to 
I'J'I 2 . For a solution with a static superfluid density, I s 
is the same everywhere along the wire, so, at a nonzero 
I s , \& cannot vanish at any point. Since vanishing of \f r 
is necessary for the phase to unwind, it is by no means 
obvious that the LAMH saddle point actually mediates 
transitions between different winding numbers. One may 
equally well imagine that the system starts at a winding 
number W, reaches the LAMH saddle point, and then 
rolls down along the negative mode, but remains all the 
time in the sector with the same winding number. 

Therefore, one is justified in looking for fluctuations 
(with time- dependent densities) that have a zero of W 
and thus allow the phase to unwind. Such fluctuations 
have been observed in numerical studies of superfluids 



flowing past obstacles |jfl, see also Ref. |J- The studies 
of Ref. Q have used boundary conditions correspond- 
ing to a fixed value of superfluid velocity away from the 
obstacle. Here we want to consider transitions between 
states with different values of superfluid velocity and so 
use periodic boundary conditions (the ring geometry) . In 
this case, the winding number in the vicinity of a uniform 
supercurrent state must be an integer. 

In a nearly uniform Bose fluid, a natural candidate for 
the critical fluctuation is one of the many time-dependent 
solutions (solitons 0) of the perfectly uniform problem. 
According to the preceding discussion, we need a field for 
which the order parameter is completely extinguished at 
some point. For a large enough length L of the ring, a 
solution having this property is 



V(x,t) = ^le-^+^^tanh: 



2& 



(1) 



Here and [i are suitably chosen constants (ji is real, 
is complex), V\ — (2n/mL)(W+ 5), W is the winding 
number of a nearby uniform state (an integer), m is the 
mass of the particle, and £i is the "healing" length, £f 2 = 
4gm| , I , i| 2 , g is the coupling constant. 

The field (JTJ solves the time-dependent Gross- 
Pitaevskii (GP) equation at zero external potential, and 
represents a depletion of density moving at velocity v±. 
The classical energy of this state is computed in the Ap- 
pendix. The momentum of this state is mv\N , where 
N is the total number of particles. Relative to the 
uniform state with winding W and superfluid velocity 
vq = 2i:W/mL, the momentum is 



AP = mv\N — mvoN = nn , 



(2) 



where n = N/ L is the linear density. In the limit L — -> oo 
with n and W fixed, the field JIJ becomes the classi- 
cal counterpart of the extremal point of the Lieb-Liniger 
spectrum [(j of the corresponding quantum problem — a 
Bose gas with a delta-function repulsion. 

Note that, while in the limit of small v\ the solution 
JU is in a certain sense close to the corresponding LAMH 
saddle point, for larger v% the two solutions are substan- 
tially different. 

We want to identify as a critical droplet mediat- 
ing transitions between supercurrent states with winding 
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numbers W and W + 1, i.e., superfluid velocities vq and 
i> 2 = 2ir(W + l)/mL. This identification shows imme- 
diately that in a nearly uniform superfluid TAPS may 
be strongly suppressed. Indeed, if the droplet is to de- 
cay into a state with winding W, the momentum AP 
needs to be transferred to another excitation branch, i.e., 
phonons, or to some external system. A look at the Lieb- 
Liniger spectrum tells us that any phonon state with mo- 
mentum AP has energy larger than the energy of the 
droplet state. 

The nearly uniform limit can be realized in a variety 
of physical systems. An interesting and perhaps some- 
what unexpected realization is a thin superconducting 
wire. In this case, the analog of Bogoliubov's phonon is 
the gapless plasmon mode describing fluctuations of 
the superconducting density. The large-scale ( "hydrody- 
namic" ) effective theory of these plasmons is a GP theory, 
with the coupling constant g — 4e 2 /C, where C is the ca- 
pacitance of the wire per unit length (see, e.g., Ref. 
This description of superconductors is specific to one di- 
mension, and the reason why it holds is that in ID the 
screening of charges is weak, so the large-scale dynamics 
is dominated by the charging energy, rather than the en- 
ergy of condensation (as in the GL theory). The order 
parameter ^ of this description is proportional (but not 
necessarily equal) to the GL order parameter. It can be 
interpreted as the field of Cooper pairs, and the entire 
description applies only at length scales larger than the 
"size" of a pair, i.e., the GL coherence length £gl- Since 
in practice £gl is much larger than the "healing" length 
£ obtained from the GP description, the size of a criti- 
cal droplet will now be determined by £gl, rather than 
£. Nevertheless, as we discuss further in the concluding 
section, it is possible to adapt some of our results to this 
case. In particular, we expect that at scales of order £gl 
disorder, typically present in superconducting wires, self- 
averages, so that as an initial approximation we can use 
the uniform limit with suitably renormalized parameters. 
At the next level of accuracy, we will need to include mo- 
mentum transfer to disorder via various mechanisms. 

Another possible realization of the nearly uniform limit 
is atomic superfluids — specifically, trapped Bose gases — 
which can in principle be prepared is such a way that 
there are no significant sources of disorder on scales 
shorter than the healing length. In this case, the above 
considerations suggest that the requirement of momen- 
tum transfer can present a major bottleneck for TAPS. 

We want to stress that, as shown in Refs. 0, ||, in 
quantum theory phase slips occur even in a perfectly 
uniform superfluid, by thermally-assisted tunneling with 
momentum transfer to phonons. The question we ad- 
dress in the present paper is whether phase slips are pos- 
sible by a classical mechanism, such as thermal activa- 
tion. This question is particularly relevant at temper- 
atures T > gn {g is the coupling constant), where the 
approximate methods used in Refs. H,@ do not apply. 

Here, we attempt to simulate TAPS numerically, us- 
ing classical equations of motion. Our simulations are 



microcanonical, i.e., there is no coupling to any exter- 
nal heat bath, and are similar in spirit to simulations of 
sphaleron transitions in the (l+l)-dimensional Abclian 
Higgs model |lCj. In our case, however, because of the 
exact solvability of the theory at zero potential , we 
take special measures to ensure appropriate population 
of the phase space (see below). 

Our main result is the complete absence, within the 
classical mechanics, of TAPS in the uniform Bose gas. 
This leaves quantum tunneling as the only mechanism of 
supercurrent decay in this case. 

We have found that TAPS appear in the presence of 
an external potential V(x), but at a rate that depends 
sensitively on the magnitude and scale of variation of V . 
When V is smooth and not too large, the rate is much 
smaller than that derived from the LAMH theory. Fluc- 
tuations that mediate TAPS in such a smooth potential 
are similar to the field JTJ. 

The paper is organized as follows. Details of the nu- 
merical procedure are described in Sects. [H] (initial con- 
ditions) and lllll (the evolution algorithm) . Numerical 
results are presented in Sect. IIVI Our conclusions are 
summarized in Sect. [Vj where we also discuss some ap- 
plications of our results. 

II. INITIAL CONDITIONS 

For the evolution equation, we have used the GP equa- 
tion 

i^ = -^d 2 x y + g\y\ 2 y + V(x)y (3) 

with periodic boundary conditions 

<S>(x + L) = V(x) . (4) 

Here \I/(cc, t) is the complex order parameter, m is the 
mass of the particles, g > is the coupling constant, 
and V(x) is an external potential. For the uniform case, 
V(x) = 0. We have set % = 1. 

Important parameters that we will often use in what 
follows are the zero-temperature speed of Bogoliubov's 
phonons 

Co = {gn mc /m) 1/2 (5) 

and the zero-temperature "healing" length 

£= (4 ff mn avc )- 1 / 2 ; (6) 

n avc is the average density of the gas. In a classical sim- 
ulation, ?i avo is obtained as the average of the density 
n = \]/t\|) over th e entire lattice: 

™avc = J J dx&V . (7) 

A precise characterization of temperature was not a 
goal of the present work. This allowed us to initialize the 
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field directly from the phonon power spectrum. Namely, 
suppose we separate the initial field W into the modulus 
(square root of the density) and the phase: 



■9{x,ti) = ^/n~[xje i6 ^ 
For the phase, we write 

2%WiX . 



9{x) 



0i(x) , 



(8) 



(9) 



where Wi is the initial winding number, 9o is an unin- 
teresting constant, which will be set to zero, and &\ is a 
fluctuation. For the density, we write 



n{x) = n avc + Sn(x) , 



(10) 



where n avc is the uniform average density, which we spec- 
ify, and Sn is a fluctuation. 

Next, the fluctuations 9\ and Sn are expressed using 
Bogoliubov's transformation 



0i(x) 
Sn(x) 



1 \ - V Zk 




ikx 



where 



u, k = co^ + fcY] 172 
Z k = g(l + k 2 e) ■ 



(11) 
(12) 



(13) 
(14) 



In our classical simulation, we take b k and bt to be ran- 
dom classical variables, whose moduli are given by the 
Bose distribution 



\b k \ 2 



(e 



n k /r 



I)- 1 



where 



O fe = ui k + vk 



(15) 



(16) 



is the quasiparticle dispersion law in the presence of a su- 
perfluid velocity v. The phase of b k is a (pseudo) random 
number uniformly distributed between and 2ir. 

In numerical work, instead of ijHJl we use the linearized 
version 



# = W^Aw 1 + 



Sn 



2n a 



(17) 



where Af is a normalization coefficient enforcing the con- 
dition 0. 

The use of the linearized Eq. (|17|l is consistent with 
the assumption Sn -C n avc of Bogoliubov's theory and 
with using the ideal-gas distribution (|15fl for phonons. 
But how far is the result from the true equilibrium of the 
interacting system and how much should we worry about 
the difference? 



To begin with, notice that for a classical field the 
notion of a "true" equilibrium applies only to the low- 
frequency component, Q k T, which is more or less clas- 
sical. For a TAPS, we expect k ~ l/£, i.e., Q k ~ gn avc , 
so that at T > gn avc a TAPS is essentially a classical 
fluctuation. For high-frequency modes, the classical ap- 
proximation is not parametrically justified, but under the 
above condition on the temperature these only provide a 
"heat bath" , and we do not expect the results to depend 
very sensitively on the model we choose for them. 

By the same token, it seems extremely unlikely that 
modest deviations from thermal equilibrium can jeopar- 
dize the qualitative results of this work (such as vanishing 
of the TAPS rate in the uniform system). 

One may, however, raise the objection that, in the uni- 
form case or for a weak potential, our way of generat- 
ing the initial conditions may not populate adequately 
the Lieb-Liniger branch, to which for instance the state 
JU belongs. To counter this objection, we subject the 
field, after it has been initialized as described above, to 
a "premixing" stage, wherein it evolves for a while in a 
relatively strong potential. After that, the potential is 
switched to the desired strength, the winding number is 
reset to the original value, and the actual evolution be- 
gins. Premixing causes further deviations from thermal 
equilibrium but these are again deemed harmless for the 
reason stated above. 

Two further comments are in order: 

(i) Bogoliubov's theory |T^ . when formulated in terms 
of density and phase, as in Ijlljl and 1)12(1. requires small- 
ness of Sn but does not require either smallness of 9± or 
a non-vanishing expectation of the order parameter. It 
is therefore applicable even in ID, cf. Ref. [Tflj . 

(ii) Within that theory, the properly subtracted vari- 
ance of the density is 



Co 



(Sn ) T - (Sn ) T=0 = - — 



dk\k\ 



1 



^9J-oo y/l + k^ eQk/T - 1 ' 

(18) 

Of main interest to us here is the temperature region 
T > gn avc . For an estimate of Sn, we will assume the 
limit 



Cf) 



2gn a 



(19) 



In this limit, the integral in l|18[) is saturated at small k: 
k ~ ^ 1 ( 1 — w 2 /cq) 1 / 2 (this infrared sensitivity is specific 
to one dimension). At these k, the Bose distribution 
can be replaced by the classical T/fifc, and upon that 
replacement the integral can be easily evaluated: 



(Sn 2 ) T 



(Sn 2 ) T=0 



(20) 



We conclude that (unless v is very close to cq) the small- 
ness of a typical density fluctuation, 5n/n avc <C 1, re- 
quires that 



T < c n a 



(21) 
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This condition is compatible with JT5J, provided the di- 
mensionless coupling g/co is chosen small. 



III. EVOLUTION ALGORITHM 

In the code, it is convenient to work with rescaled, 
dimensionless variables, since this reduces the number 
of parameters that need to be input. A natural unit of 
length is the healing length £, and of time — the ratio 



1 



(22) 



Co 2gn avc 

So, we rescale variables in the GP equation as follows: 

xli , (23) 

2gn avc t , (24) 

V/2gn avo ■ (25) 



.r 
t 

V 



In this section, and in the plots, we use the same letter 
for a rescaled quantity as in the rest of the paper for the 
original one. We also rescale the field \& into 



(26) 



so that the average rescaled density is equal to 1: 
(1/L) J i/jtydx = 1. 

In the rescaled variables, the GP equation becomes 



dip 

'at 



(27) 



For the purpose of numerical evolution, we associate to 
the right-hand side of l|27|l an operator U that updates 
the field ^p n {x) at the n-th time step: 



r/i n+ i(x) = U[^ n (x),At] ; 



(28) 



At = f n +i — t n . For this work, we have used an operator- 
splitting algorithm, in which the update is made in three 
steps: 

Uty(x),At]=U 1 [U 2 [U 1 ty(x),±At},At],^At] . (29) 
where U\ corresponds to the equation 



i^ = ~\i>\ 2 iP + V(x)iP 



and U2 to the equation 



"at = 



(30) 



(31) 



To solve these equations, the field was discretized on 
a spatial lattice with a uniform step Ax and periodic 
boundary conditions. Eq. I|3UI) is local in space, so U± 
can be defined site-by-site: 



M-l 

U x [i>{x),At] = J] 
3=0 



(32) 



We have used 



1 + (To 



(33) 



with Uj — ^[^\tp\ 2 + V(xj)]At. The operator C/ 2 is non- 
local in x-space but local in fc-space; we have used 



M-l 



U 2 [i/>(x),At] = J] v k (^,At) 



(34) 



fe=0 



where ip F is the fast Fourier transform (FFT) of %p, and 



v k (^ F ,At) 



1 - Pk 1 F 



l + pk 



(35) 



with pk = iifc 2 . So, there are two FFTs (one direct and 
one inverse) at each time step. 

The resulting algorithm has second-order accuracy in 
time. It conserves the number of particles exactly, while 
the energy non-conservation is controlled by the time- 
step At. 



IV. NUMERICAL RESULTS 

We present results from simulations on an M — 1024 
lattice with the following values of the dimensionless pa- 
rameters: L/£ = 79, T/2gn avc — 1.2, and g/co = 0.1 (a 
weak coupling). We have used the periodic potential 



V(x) = y o sin(27rx/0 



(36) 



For periodicity, me must have I = L/q, where q is an 
integer. We have used q = 15, so that //£ w 5.3. "Pre- 
mixing" was done with Vo/2gn llvc — 0.3 and lasted for 
At = 200 (in our dimensionless units). The premixing 
stage is not shown in the plots, i.e., t = corresponds to 
the beginning of the actual evolution with a smaller po- 
tential. The winding number changes during premixing, 
but at t = it is reset to the original value Wi = 3. 

In Fig. ^ we show the evolution of the winding num- 
ber with and without a potential. In the uniform case 
(Vo = 0), the winding number fluctuates but always 
quickly returns back to the initial value W = 3. There is 
no overall relaxation. This behavior persisted as we went 
to larger values of the temperature. Based on these and 
similar other results, we conclude that, in the uniform 
ID superfluid, supercurrent is classically stable for any 
value of superfluid velocity satisfying Landau's criterion. 
Thus, the only mechanism of supercurrent decay in this 
case is a quantum effect — the phonon-assisted tunneling 
considered in Re fs. 

In the presence of a (weak) potential, while there are 
still many unsuccessful phase-slip attempts, the winding 
number eventually relaxes towards zero. 

These results may be compared to those of numer- 
ical simulations of sphaleron transition in the (1+1)- 
dimensional Abelian Higgs model ^(|- Sphalerons are 
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FIG. 1: Evolution of the winding number for Vo/2gn { 
0.03 and ^0 = 0. 
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FIG. 2: Power spectrum of the field ip for k > in the uniform 
system at the beginning and near the end of the simulation. 



very similar to TAPS except that they connect states 
with zero current, so there is no issue of momentum trans- 
fer. In the simulations of Ref. ^(|> they readily occur in 
the absence of any external potential. 

In Fig. |3 we show the power spectrum of ip in the case 
V{) = at the beginning and near the end of the simu- 
lation (only positive k are shown; results for k < are 
similar). We see that the spectrum exhibits remarkable 
stability. Although one expects that eventually fluctu- 
ations will propagate to the ultraviolet (a manifestation 
of the Ray leigh- Jeans problem of classical statistics), this 
clearly does not happen on the timescale of our simula- 
tion. 

In Fig. |3 we plot the density profiles of the fluctuation 
mediating one of the phase slips in Fig. for Vo 7^ 0, 
at three different moments of time. Comparing these 
profiles with Eq. 0J, we find that the spatial size of the 
droplet is about 2.5 times smaller than predicted by that 
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FIG. 3: The "flying droplet": a moving depletion of density 
that mediates a phase slip in the presence of a potential. 

equation. This is not surprising given that, in the present 
simulation, fluctuations of density are relatively large: 
Sn is about half of n avc . Since the typical wavelength 
of these fluctuations is of the same order as the size of 
the droplet [cf. comment (ii) in Sect. at Sn ~ rt avc 
we loose the notion of an effective long-range theory to 
which could be a solution. 

This interpretation is confirmed by going to smaller T 
(and stronger potentials), for which Sn/n we is smaller, 
and the droplet size is expected to be closer to the value 
predicted by Eq. (Q. For example, for T/2gn avc = 0.9 
and Vb/2gn avc =0.1, the difference reduces to a factor of 
about 1.6. 



V. DISCUSSION 

The main result of this paper is that in a uniform 
ID Bose gas classical phase slips are completely blocked 
out. We find this result nontrivial and even surprising, 
given that momentum conservation, while characteristic 
of the uniform system, by itself does not prohibit phase 
slips: momentum released from the supercurrent can be 
absorbed by phonons. Indeed, quantum-mechanically, 
phase slips are possible even in the perfectly uniform sys- 
tem Hit What we have shown here, then, is that quan- 
tum effects remain the only source of supercurrent decay 
in the uniform case. 

Our results, both for the uniform system and in the 
presence of a potential, are consistent with identifying the 
moving droplet Q as the fluctuation mediating thermal 
phase slips. As we have already noted, at nonzero v% it 
is distinct from the LAMH saddle point. 

The nearly uniform limit can presumably be realized 
in trapped Bose gases. Another system where, as already 
mentioned in the Introduction, the droplet Q may play a 
role is a thin superconducting wire. In this case, the field 
$ describes a fluid of Cooper pairs, and for the reasons 
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already indicated — the weak screening in ID and the re- 
sulting special role of the charging energy — the GP equa- 
tion is well justified for description of long-wavelength 
plasma oscillations — the gapless plasmon mode . How- 
ever, this description breaks down at scales of order of 
the "size" of a pair, i.e., the Ginzburg-Landau coherence 
length £gl- Since in practice £gl is much larger than 
the healing length £ obtained from the GP description, 
the size of the critical droplet will now be determined by 
£gl> rather than £. In the following discussion, we will 
need only the energy of droplet, which can be calculated 
using the Ginzburg-Landau free energy. The latter dif- 
fers from the GP Hamiltonian I|A1(I only by the value of 
the quartic coupling g and by the presence of a 1^1 2 term, 
which in the GP case would correspond to a chemical po- 
tential. Although the activation energies l|A7(l . (|A8|) are 
computed in the Appendix at a fixed number of parti- 
cles, the same expressions also apply at a fixed chemical 
potential. As for the coupling g, it will be sufficient to 
consider it as a phenomenological parameter. This allows 
us to apply the expressions obtained in the Appendix to 
the case of a superconducting wire. 

Since a moving droplet can now transfer momentum 
to normal electrons and, through them, to the disorder 
potential, we expect that for sufficiently strong disorder 
the TAPS rate is determined simply by the droplet's acti- 
vation energy. The activation energies for transitions de- 
creasing and increasing the winding number by one unit 
are given by Eqs. I|A7|I . I|A8[1 . Unfortunately, the form 
of the I—V curve obtained from these expressions differs 
from that obtained in the LAMH theory only for currents 
comparable to the GL critical current. This makes it dif- 
ficult to draw an experimental distinction between the 
activation energies (|A7|I . (|A8|) and their LAMH counter- 
parts. If, at smaller currents, we regard the last term in 
Eqs. (|A7|I . I|A8I) as current-independent, the voltage drop 
found from these equations is proportional to sinh(/ s //o) 
where I s is the supercurrent, and Iq = 2eT /it. This form 
of the nonlinear I—V curve, the same as in the LAMH 
theory in the equivalent limit 0] , has been recently found 
to be in good agreement with experiment |l4j| . 

Finally, we mention one technical result of our work, 
namely, the extent to which it turned out possible to 
simulate thermal field theory classically. We refer here 
to the remarkable stability of the classical power spec- 
trum against spreading towards the ultraviolet, cf. Fig. 
|2| This is similar to behavior observed in three dimen- 
sions, in simulations of highly nonthermal states pro- 
duced by parametric amplification [15J . One may wonder 
if this property holds more generally, so that the classical 
model of a thermal state we have used here can be ap- 
plied to other (weakly-coupled) systems where the high- 
frequency modes provide an internal "heat bath" for the 
low- frequency, classical component. 

After this work was completed, we have learned of a 
paper by Polkovnikov et al. |l6j| who studied numerically 
classical decay of supercurrent in an optical lattice. Their 



calculation corresponds to a strong periodic potential — 
the limit opposite to ours. 

The author thanks V. Ambegaokar, L. Pitaevskii, A. 
Rogachev, and A. Zubarev for useful discussions. This 
work was supported in part by the U.S. Department of 
Energy through Grant DE-FG02-91ER40681 (Task B). 



APPENDIX A: ENERGY OF A MOVING 
DROPLET 

Here we compute the energy of the field Q), rela- 
tive to the energy of the neighboring uniform states: 
one with winding number W and supcrfluid velocity 
vq = 2TrW/mL, and the other with winding W + 1 and 
superfluid velocity V2 — 2n(W + l)/mL. 

We start with the Hamiltonian corresponding to Eq. 
© with V = 0, 



dx [ — \d x V 
2m 



2 + !w 4 



(Al) 



and compute the energy at a fixed total number of par- 
ticles N. So, there is no need to include a chemical po- 
tential. 

Substituting Q into I|A1(I and integrating over x, we 
obtain (to exponential accuracy in £i/L) 

£i = y^ + f|*i| 4 (L-^i), (A2) 

where £i is the healing length defined after Eq. Q . The 
particle number is 



N 



J dx^V = (£-4fi)|#i| 



(A3) 



The energies of the neighboring uniform states are 



£o,2 = y< 2 iV+§|*o| 4 £ 



(A4) 



where |<F | 2 = N/L. To the leading order in 1/L, £i s» £, 



2 2 
"1 ~ v 2 



2ttvi 
mL 

l*i| 2 -|*o| 2 w ^|*o| 2 



(we allow for the possibility that W oc L). Then, 



(A5) 
(A6) 



S 1 -S «7r Ul |^ | 2 + - 3 el*o| 4 , (A7) 

Ei--E 2 «-7roi|* | 2 + |<tf|*o| 4 . (A8) 

These expressions coincide with the LAMH activation 
energies |2j expanded to the first order in supercurrent. 
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